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ABSTRACT 

We show that the mean values and distributions of the time ending the shallow 
decay of the light curve of the X-ray afterglow of long gamma ray bursts (GRBs), 
the equivalent isotropic energy in the X-ray afterglow up to that time and the 
equivalent isotropic GRB energy, as well as the correlations between them, are 
precisely those predicted by the cannonball (CB) model of GRBs. Correlations 
between prompt and afterglow observables are important in that they test the 
overall consistency of a GRB model. In the CB model, the prompt and afterglow 
spectra, the endtime, the complex canonical shape of the X-ray afterglows and 
the correlations between GRB observables are not surprises, but predictions. 



1. Introduction 

In a relatively brief time, the observations made or triggered by the Burst Alert Telescope 
(BAT), the X-Ray Telescope (XRT), and the UVOR telescope, aboard the Swift satellite, 
have gathered a wealth of new information on Gamma Ray Bursts (GRBs), in particular on 
the early X-ray and optical afterglows (AGs) of long-duration GRBs. These observations pose 
severe problems to the generally accepted 'Fireball' model of GRBs (see, e.g. Meszaros 2006), 
whose 'microphysics' (see, e.g. Panaitescu et al. 2006), reliance on shocks (see, e.g. Kumar 
et al. 2007), and correlations based on the 'jet-opening angle' (see, e.g. Sato et al. 2007; 
Burrows & Racusin 2007), may have to be abandoned. 

The said recent observations agree remarkably well with the predictions of the 'Cannon 
Ball' (CB) model (Dar & De Riijula 2004; Dado, Dar & De Rujula 2002; Dado, Dar & De 
Rujula 2003, hereafter DD04; DDD02; DDD03, respectively). Some examples are given in 
Fig. [H The predicted lightcurve of the X-ray AG afterglow is shown in Fig. [T^l for the fireball 
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(see e.g. Maiorano et al. 2005) and CB (DDD02) models. For GRBs with an approximately 
constant circumburst density distribution, the well-observed X-ray AGs have a 'canonical 
behaviour' (e.g. Nousek et al. 2006; Zhang et al. 2006; O'brein et al. 2006), in impressive 
agreement with the CB-model predictions, as in the example of Fig. [lb. The evolution of 
the AG around the time, T a , ending the 'shallow phase' is predicted to be achromatic in the 
optical to X-ray range (DDD02), as can be seen in Fig. [lb and its comparison with Fig. [Tbl. 
This 'achromaticity' does not extend to the radio domain (DDD03). Many correlations 
between (prompt) GRB observables have been studied with the help of the new data on 
GRBs of measured red-shift z (e.g. Schaefer 2006). All of the successful correlations are 
simple predictions of the CB model. One example, perhaps the best known, is the correlation 
between the spectral 'peak' energy and the total (bolometric) isotropic-equivalent energy, 
shown in Fig. [Tbl along with our prediction (Dado at al. 2007 and references therein). 

Willinger et al. (2007) have studied and tabulated recent data on T a , the time ending 
the shallow decay of the X-ray light curves of long GRBs. In a paper the title of whose first 
version was the same as ours, Nava et al. (2007) have investigated the correlation between 
T' a = T a / (1 + z), the endtime in the source's rest frame, and the prompt GRB energy, as well 
as the correlation between T' a and the energy in the X-ray plateau phase, integrated up to 
T a . The values and distributions of these three quantities, as well as their correlations, are 
important in ascertaining the global validity of a GRB model, among other things because 
they test the consistency of the description of the prompt and afterglow phases. In this paper 
we show that the observations are in precise agreement with the CB-model's predictions. To 
do so, we gather, in Sections [31 SI and predictions from various of our papers, in a manner 
and order adequate to the discussion of the subject at hand. The predictions are derived for 
'typical' or average values of the parameters, all chosen as in our earlier work, which referred 
mainly to pre-Swift observations. The incidence and explicit origin of the variability around 
the typical cases is also discussed. The results are summarized in Fig. [H which demonstrates 
that the central expectation, variability and correlations are all as predicted. 



2. The CB model 

In the CB model (Dar & De Riijula 2000, DD2004; DDD02, DDD03), long-duration 
GRBs and their AGs are produced by bipolar jets of CBs, ejected in core-collapse SN ex- 
plosions (Dar & Plaga 1999). An accretion disk is hypothesized to be produced around 
the newly formed compact object, either by stellar material originally close to the surface 
of the imploding core and left behind by the explosion-generating outgoing shock, or by 
more distant stellar matter falling back after its passage (De Riijula 1987). As observed 
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in micro quasars, each time part of the disk falls abruptly onto the compact object, a pair 
of CBs made of ordinary plasma are emitted with high bulk-motion Lorentz factors, 7, in 
opposite directions along the rotation axis, wherefrom matter has already fallen onto the 
compact object, due to lack of rotational support. The 7-rays of a single pulse in a GRB 
are produced as a CB coasts through the SN glory -the SN light scattered by the SN and 
pre-SN ejecta. The electrons enclosed in the CB Compton up-scatter (Shaviv & Dar 1995) 
glory's photons to GRB energies. 

Each pulse of a GRB corresponds to one CB. The emission times of the individual 
CBs reflect the chaotic accretion process and are not predictable. At the moment, neither 
are the characteristic baryon number and Lorentz factor of CBs, which can be inferred 
from the analysis of GRB afterglows (DDD02; DDD03; DD04). Given this information, 
two other 'priors' (the typical early luminosity of a core-collapse supernova and the typical 
density distribution of the parent star's wind- fed circumburst material), and a single extra 
hypothesis (that the wind's column density in the 'polar' directions is significantly smaller 
than average) all observed properties of the GRB pulses can be successfully predicted without 
the introduction of any ad-hoc parameters (Dar & De Riijula 2004, thereafter DD04). 

The spectral energy density, F u (t), of the X-ray emission of a GRB has two phases. The 
first is very rapidly declining X-ray emission and dominated by the late-time tail of the GRB 
pulses (DD04) and/or by line emission and thermal bremsstrahlung from the CBs (DDD02; 
Dado et al. 2006). In a second phase, synchrotron radiation from swept-in ISM electrons 
spiraling in the CBs' enclosed magnetic field takes over. This second phase has a 'plateau': 
a shallow time-dependence lasting until the CBs decelerate significantly in their collisions 
with the interstellar medium, after which F v {t) bends into an asymptotic power-law decline 
~ t~ 16 . On this basis we were able to predict (DDD02) the 'canonical' behaviour of X-ray 
AGs, observed by Swift (Dado et al. 2006). 

In the CB model, three times characterize the evolution of F v (i). The first is the time 
at which synchrotron radiation begins to dominate. The second is the time at which the 
injection bend of the electron energy spectrum within a CB crosses a particular frequency, 
and corresponds to a strongly chromatic change in F v (t) (DDD02; DDD03). The third time, 
the typical deceleration time, corresponds to an achromatic change in F v {t) (DDD02) and 
is the subject of this paper. 
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The deceleration time 



Let 8 = 0(1 mrad) be the typical viewing angle of an observer of a CB that moves with 
a typical Lorentz factor 7 = 0(1O 3 ). Let 5 = O(10 3 ) be the corresponding Doppler factor: 

1 2 7 



7 (l-/3cos0) l + 7 2 # 2 
where the approximation is excellent for ^ C 1 and 7 ^> 1. 



(1) 



A CB is assumed to contain a tangled magnetic field in equipartition with the ISM 
protons that enter it. As it ploughs through the ionized ISM, a CB gathers and scatters its 
constituent protons. The re-emitted protons exert an inward pressure on the CB, countering 



its expansion. Let n 



n.„ 



n e be the number density in a dominantly hydrogenic ISM. In 



the approximation of isotropic re-emission in the CB's rest frame and a constant n, one finds 
that within minutes of observer's time t, a CB reaches a nearly-constant 'coasting' radius R. 
Subsequently, j(t) obeys: 
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with k = 1(0) depending on whether the re-emitted ISM particles are a small (large) fraction 
of the intercepted ones. This dichotomy is too small to detect in the study of AGs, but the 
case k — 1, which we adopt, is favoured by the CB model of Cosmic Rays (Dar & De Rujula, 
2006). 

To specify the CB-model's prediction for an end-time, T' a , paraphrasing the one defined 
by Willinger et al. (2007) and Nava et al. (2007), let us derive the time at which the X-ray AG 
is smaller by a factor of two than the extrapolation from its previous shallow behaviour, and 
let us refer to the typical parameters of observed GRBs, for which 7o#~l, and 5(t) ~7(t) 
in the shallow phase (DDD02). The X-ray AG, as we shall recall in Section [51 Eq. ( TT~2T) . 
behaves as F x (t) (x r y(t) 6A , so that we are demanding that [7(t)/7 ] 6 ' 4 ~l/2. Insert this into 
the left hand side of Eq. (T5]), with k = 1 and #70 = 1, to conclude that the typical end-time is 
t = 1.026 to- Nava et al. (2007) correct this time for the cosmological redshift; according to 
Eq. (J2), its predicted value is: 
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where we have normalized to typical CB-model parameters and the 'variability' around them 
is governed by the combination of parameters Vr (DDD03; Dar & De Rujula 2006). 
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4. The isotropic energy in the prompt GRB 

In the CB model the isotropic (or spherical equivalent) energy, E™°, of a GRB, is (DD04): 

- - ~ ~ ( 2 - 8 x 1053 ef g) v *> 

' be V 47T 

V = U) 
E ~ 10 9 1% 4.5 10 50 ' 1 ] 

where L SN is the mean SN optical luminosity just prior to the ejection of CBs, N CB is 
the number of CBs in the jet, iV& is their mean baryon number, (3 S is the comoving early 
expansion velocity of a CB (in units of c/v3), and a T is the Thomson cross section. The 
early SN luminosity required to produce the mean isotropic energy, ~ 4x 10 53 erg, of 
ordinary long GRBs is L|?™ ~5x 10 42 ergs _1 , the estimated early luminosity of SN1998bw. 
All quantities in Eq. are normalized to their typical CB-model values. For N CB we took 
the result of a recent careful analysis of the number of significant peaks in a GRB light curve 
(Schaefer 2006) rather than the one we previously adopted (N CB = 6, DD04). 

Nava et al. (2007) choose to present their results in terms of the prompt isotropic energy 
in the 15-150 keV domain. To restrict the 'bolometric' result of Eq. (j3J) to a fixed-energy 
bracket, we must recall the prediction of our model for a GRB's spectral shape (DD04). The 
photons of the glory's light that a GRB Compton-upscatters have a thin-bremsstrahlung 
spectrum dN/dEi oc (1/ Ei) a Exp[—Ei/Ti], with a ~ 1 and Tj ~ 1 eV. The bulk of these 
electrons are comoving with the CB and Lorentz- and Doppler- boost the target light to a 
spectrum of the same shape, and 'temperature': 

(l + z)T~ -Ti-f 6 (1 + cos 9i) ~ 1.3MeV(l + cos^), (5) 
3 

where #j is the angle of incidence of a glory's photon into the CB, in the SN rest system. 
A very tiny fraction of the moving electrons is due to 'knock-on', or is 'Fermi- accelerated' 
within the CB, in both cases to a spectrum (in the CB's rest frame) dN e /dE e oc E~ p , with 
p ~ 2.2. The complete prompt 7 spectral distribution, upscattered by the CB's comoving 
and knock-on electrons (DD04), is: 

^-(i)"^ /r+6(1 -^ ,T) (i)" 2 (6) 

For b = 0(1), and a and p in their expected range, the above spectrum is uncannily similar 
to the phenomenological 'Band' spectrum (Band et al. 1993). The 'peak energy' of the 
prompt spectrum is: 

(i+-)^-4^" (5ookev) i^i^v- (7) 
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The average redshift of Swift GRBs, and of the ones discussed here, is (z) =2.8. The 
fraction / 7 of the bolometric lying in the 15-150 keV range is the ratio of J E dN^ in the 
range (1 + z) x (15-150) keV, to the same integral from to oo: / 7 ~ 0.106 for all parameters 
at their central values and (l + cos^) = 1/2 (a semitransparent glory, DD04). Our prediction 
is then: 

£f [15-150 keV] = / 7 £f = (2.9xl0 52 erg) V E , (8) 

with Ve as in Eq. (j4j). As the rest of our results, / 7 is computed for 'typical parameters', 
corresponding to a relatively large E p value and the concomitant large bolometric corrections. 
For many post-Swift GRBs the bolometric correction would be smaller. 



5. The isotropic energy in the X-ray plateau phase 

In the plateau phase and thereafter, the CB-model's AG is due to synchrotron emission 
by the electrons continuously entering a CB from the interstellar medium (ISM) it sweeps. 
Above observer's radio frequencies, and in the CB's rest system, the synchrotron radiation 
has a (normalized) spectral shape (DDD03): 
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where the 'injection bend ' frequency corresponds to the energy, Ef, = 
ISM electrons enter the CB at the time when its Lorentz factor is 7(i). 



frequency v h 



in the CB's (observer's) frame is: 
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The typical frequency in the parenthesis is equivalent to an energy of 3.9 eV. This is always 



below the X-ray domain, so that the corresponding X-ray spectrum has a 



,-i.i 



shape. 



But, occasionally, at times of order 1 day or less, the observed optical frequencies are above 



obs 



so that the optical spectrum varies from 



~z/-°- 5 to 



•v 



-i.i 



, producing a chromatic break 



occurring in the optical AG but not in the X-ray one (DDD03). 

In a CB's rest frame, the energy flux density in the optical to X-ray domain is: 



^CB K t] — V 71 R nm e C j(t) /sync(M), 



IT 
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where rj is the fraction of ISM electrons that enter the CB and radiate there the bulk of their 
incident energy, and / sync is as in Eq. (jUJ). The AG spectral energy density F obs seen by a 
cosmological observer at a redshift z, is: 



F obe [v,t] ~N C 



CB 



[l + z)u 6{t)t 



l + z 



5(t) 



cx 7 (t) 2 - 3 5(t) 
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(12) 



where F CB is as in Eq. ffTTj) . and Dl is the luminosity distance. As announced in the derivation 
of Eq. (j3J), F obs oc 7 6 ' 4 in the shallow phase of the X-ray afterglow. 

With use of the spectrum of Eqs. (l9"fTU|) we can define a fraction / of the spectral energy 
in the 15-150 keV X-ray band. For p = 2.2: 
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Gathering the above results, we obtain for the equivalent isotropic energy per unit observer 
time in the specified X-ray range: 



[7(f)] 


2.3 


\s(ty 


4.1 


3.8 


i.i 


n 


1.05 


R 


2 




. 103 . 




_10 3 _ 




[(l + z)_ 




.10~ 2 cm -3 . 




_10 14 cm_ 




. 4 - 5 . 



F x [15-150keV] = f F x = N CB rj f n R nm e c r fbj (1 + z) ~ (9.7 x 10 47 ergs 
V F (t) 

where we used (1 + z) — 3.8 for the average redshift of GRBs detected by Swift. 
The corresponding integrated X-ray energy in the plateau is: 

1 



X 



^F x (t = 0)T a = (6.8 x 10 50 erg)V r V>(0), 



VV F 
(14) 



(15) 



where we took rj ~ 1, Vt and Vp are as in Eqs. (jSJ) and ( JT4|) . and the factor 1/2 reflects the 
fact that (as can be seen in a plot which is not a logarithmic) most of the plateau, as defined 
here, extends in the domain where the AG light curve has a value ~ 1/2 of its initial value. 



6. Results 

6.1. Afterglow versus prompt bolometric energies 

A simple and crucial test of models of GRBs is the predicted ratio of the bolometric 
energy in a GRB's afterglow up to the end of the plateau phase (essentially all of the AG's 
energy) and the total energy in the GRB's prompt 7 rays. According to Eqs. ( fT4l) and ( |T5l) . 



- 8- 



the CB-model expectation is: 

i?[AG/GRB] = ^ J m e c 6 = 0.08 V R 




Vr ~ TTe^lJs^^^- (16) 

This ratio is rather 'clean': it establishes a link between the late and prompt emissions which 
is independent of the number of CBs, of their radii, of the density of the ISM in which they 
travel, and weakly dependent on their baryon number. It very naturally explains why the 
observed ratios are typically of the order of a few percent. 



6.2. Central and typical values 

Our main results are Eq. ([3]) for the time ending the shallow X-ray AG decay, Eq. (jSJ) 
for E™°[15 — 150 keV], a GRB's prompt isotropic-equivalent energy in the specified interval, 
and Eq. (ITS"]) for the isotropic energy, P^°, in that energy interval of the X-ray AG, integrated 
in time up to the end of the plateau. The 'typical' parameters underlying these results are 
based on the analysis of pre-Swift GRBs, and reflect the domain wherein GRBs, in the past 
and with less performing satellites, it was easiest to detect GRBs. 

The predicted central expectation for T' a is shown in both parts of Fig. [2] as a horizontal 
line. The central expectations for E^°[15 — 150 keV] and P so are the vertical lines in the 
upper and lower part of the figure, respectively. The 'sweet spot' drawn as the small ellipses 
corresponds to letting 7 and 5 vary in the observed narrow domain wherein most pre-Swift 
GRBs lied (DDD02; DD04). The larger dotted ellipses are drawn by allowing the combina- 
tions of variability parameters, V in Eqs. (pH El US]), vary by about an order of magnitude 
around the small ellipse. As expected, most Swift GRBs with relatively large and measured 
peak energy (the stars) and relatively large prompt and AG isotropic energies, are within 
the dotted ellipses. Most of the extra points (the dots) may in the past have been classified 
as XRFs (they have low 'peak energy', E p ). We discuss them in the next subsection. 

The green lines in Fig. [2] show the correlations expected for typical parameters. For 
them, the 7 and S dependences of the relevant quantities are T' a oc -y^ 3 , E™° oc 5q and 
Jf° oc 7o~ ' 7 ^ 1 , so that T' a oc l/E iso and T' a oc (ijp)" 3 / 3 - 4 for S w 70. 
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6.3. Distributions and correlations 

In the CB model, XRFs are the same as GRBs, but observed at a relatively large angle, 9 
(or a particularly small 7), implying a small 7 5 (DD04; Dado et al. 2004). Thus, XRFs have 
a relatively small spectral peak energy and a small prompt isotropic energy. The explicit 
proportionality factors in the relations E p oc 7 5 and E™° oc 5 3 are given by Eqs. (17114"]) . 
Consider them fixed at their typical values. The typical (7, 5) domain of observable GRBs 
is then the one shown in Fig. [31 The observed values of 7 are fairly narrowly distributed 
around 7~10 3 (DDD02, DD04), as in the blue strip of the figure. The (7, 5) domain is also 
limited by a minimum observable isotropic energy or fluence (both oc 5 3 ), by a minimum 
observable peak energy, and by the line 9 = or by a line corresponding to a minimum 
fixed 9, if one takes into account that phase space for observability diminishes as 9 — > 0. 
The elliptical 'sweet spot' in Fig. [3] is the region wherein GRBs are most easily detectable, 
particularly in pre-Swift times. X-ray Flashes populate the region labeled XRF in the figure, 
above the fixed 76? line or to the left of the fixed E p line. We interpret most of the dotted 
points in Figs. [3] and H] as cases for which 7 and S lie in the 'XRF domain' of Fig. [31 

The continuous red lines in Fig. [Hare the contours of the blue domain of Fig. [31 projected 
into the [E™°, T^] plane (top) and the [£>x°, T^] plane (bottom). The projectors' are the 
corresponding functions of 7 and 5, e.g. ^(7, 5), as in Eq. ([3]). The dotted red lines are drawn 
by 'moving' the red contour about its 'central' position in the planes, by approximately one 
order of magnitude, once again to reflect the variability of parameters other than 7 and 9 
(or 7 and 5). The dotted red lines satisfactorily describe the location and distribution of 
the Swift data. The green lines are the predicted trend of the correlations, which for the 
ensemble of the data (stars and points) interpolates between two power laws, as in Fig. [TH, 
and as discussed in detail for this and many other correlations in Dado et al. 2007. 

7. Conclusions 

We have analised data on two afterglow observables (the time ending the shallow decay 
of X-ray AGs and the integrated isotropic energy up to that point) as well as a prompt- 
GRB observable (the isotropic energy). To do so, we have simply reported the theoretical 
expectations of the CB model, developed in previous papers. The predictions include the 
magnitudes of these quantities, the explicit dependence on the parameters that govern their 
case-by-case variability, and the spectral shapes of the prompt and afterglow phases. 

The results can best be summarized by looking at Fig. HI whose data points are those in 
the corresponding figure in Nava et al. (2007) who discuss, among others, the same subject. 
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The predictions of the CB model, predating the Swift data, are in excellent agreement with 
the observations. The data are centered and distributed as expected. Their correlations are 
also the expected ones, though, since the data points do not span a large number of orders of 
magnitude, they are not as remarkable as for other correlations, such as the one in Fig. [Tfci. 

In our opinion, the main novelty in the paper of Nava et al. (2007) is the discussion of 
correlations between prompt-GRB and GRB-afterglow observables. These test the ensemble 
and coherence of a GRB-model's ingredients. In the CB model the prompt 7 rays are of 
Compton origin, while the AG light is dominated by synchrotron radiation. Unlike fireball 
models based on very different physics, the CB model never had an 'energy crisis' (see 
e.g. Piran 2000) in the relation between the total energies in the prompt and afterglow phases, 
or a problem with the prompt spectrum (see, e.g. Ghisellini 2003). The time ending the 
shallow afterglow phase is a 'deceleration time' of the cannon-balls, unrelated to the opening 
jet-angle of fireball models. The CB model provides very simple, predictive and successful 
descriptions of the prompt (DD04, Dado et al. 2007) and afterglow phases (DDD02; DD03; 
Dado et al. 2006; Dado et al., in preparation). We are not surprised that the model is also 
successful in the detailed description of the distributions of prompt and afterglow observables, 
and of their correlations. 

Acknowledgment: This research was supported in part by the Asher Space Research 
Fund at the Technion. 
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Fig. 1 — Top left: Pre-Swift predictions for the 2-10 keV X-ray AG in the CB (DDD02) 
and fireball (Maiorano et al. 2005) models, compared to data for GRB 990123. T a is the 
time ending the plateau phase. Top right: Broad band optical data on GRB 990123, fit in 
the CB model (DDD03). The evolution is achromatic all the way up to the X-ray energies. 
Bottom left: Comparison between the CB model prediction and the canonical 0.2-10 keV 
X-ray light curve of GRB 050315 (Vaughan et al. 2006) Bottom right: The (E p , E*°) 
correlation, compared with its predicted trend in the CB model (Dado et al. 2007). The 
crossed red lines are the predicted typical or average values, see Eqs. ( HfTj) . 
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Fig. 2. — The central and typical values, discussed in Section 16. 2\ for the time ending the 
X-ray plateau, plotted against the isotropic energies of the prompt 7 rays and of the X-ray 
AG up to that time. The data are those gathered by Nava et al. (2007). The crossing 
lines are the predictions of Eqs. (l3ll8]TT5l) . Within the small ellipse, the parameters 70 and 
So are allowed to range in the small domain in which most pre-Swift GRBs gathered. The 
larger ellipse allows for the relevant combinations of the other case-by-case parameters to 
vary by about one order of magnitude. The thick (green) lines are the expected trend of the 
correlations. 
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Fig. 3. — The domain of (5, 7) values. The green elliptical spot labeled 'GRB' is the area 
wherein most pre-Swift GRBs were observed. The region to its left has relatively small 
(large) S (6) values, corresponds to relatively small and E p , and is labelled 'XRF'. The 
blue horizontal band is limited above and below, reflecting the narrow distribution of the 7 
values of observable CBs (DDD02; DD04). 
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Fig. 4. — The time ending the X-ray plateau, plotted against the isotropic energies of the 
prompt 7 rays and of the X-ray AG up to that time. The data, crossing lines and ellipses 
are as in Fig. [2j The red line in the top figure is the projection of the contour of the blue 
(7, 5) domain of Fig. [3] onto the (T' a , E™°) plane. It encompasses the area in which GRBs 
and XRFs are expected to lie, for all parameters set to their central values, but for 7 and 
5. The red dashed contour is obtained by letting the rest of the relevant combinations of 
parameters vary by about one order of magnitude. The green thick lines show the expected 
trend of the correlations. The lower figure is built in the same way, in the (T' a , E 1 ^) plane. 



